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One of the hallmarks of spintronics is the control of magnetic moments by electric fields enabled 
by strong spin-orbit interaction (SOI) in semiconductors. A powerful way of manipulating spins in 
such structures is electric dipole induced spin resonance (EDSR), where the radio-frequency fields 
driving the spins are electric, and not magnetic like in standard paramagnetic resonance. Here, we 
present a theoretical study of EDSR for a two-dimensional electron gas in the presence of disorder 
where random impurities not only determine the electric resistance but also the spin dynamics via 
SOI. Considering a specific geometry with the electric and magnetic fields parallel and in-plane, 
we show that the magnetization develops an out-of-plane component at resonance which survives 
the presence of disorder. We also discuss the spin Hall current generated by EDSR. These results 
are derived in a diagrammatic approach with the dominant effects coming from the spin vertex 
correction, and the optimal parameter regime for observation is identified. 

PACS numbers: 73.23.-b, 73.21.Fg, 76.30.-v, 72.25.Rb 



The field of spintronicsi*^ focuses on the interplay be- 
tween spin and charge degrees of freedom of the elec- 
tron. The relativistic effects responsible for the cou- 
pling between spin and orbital motion can be strongly 
enhanced in solids due to band structure effects, with 
III-V semiconductors showing a particularly strong spin- 
orbit interaction (SOI) resulting in zero-field spin split- 
tings. For instance, bulk inversion asymmetry gives rise 
to Dresselhaus SOI)^ while structural inversion asymme- 
try occurring in heterostructures gives rise to Rashba 
SOIi^ The strength of such SOIs can be varied over 
a wide range which offers the advantage to control 
magnetic moments with electric fields. A well-known 
and particularly powerful way of manipulating spins 
in such structures is electric dipole induced spin reso- 
nance (EDSR) r'SiSiSiifliii where the radio frequency (rf) 
fields coherently driving the spins are electric, and not 
magnetic like in standard paramagnetic resonance^^ii^ 
The advent of materials with tailored SOli has sparked 
intense interest in a variety of spin orbit effects and 
its applications such as spin currentsji^iiSiiSiiLiS gate- 
controlled SOI effects, ^iSiSfliSiiSS spin relaxation in quan- 
tum dots,'^'^'^"'~ zitterbewegungf2& spin-based quan- 
tum information processing^*2i and, in particular, 
EDSR)^i^°i" which is the focus of this work. 

Experimental indication for EDSR was recently re- 
ported for semiconductor epilayers using an in-plane elec- 
tric rf field^ (with the magnetic field applied out of 
plane). In this geometry, spin coupling to the elec- 
tric field is much stronger than when the electric field 
is applied along the growth direction.-- Further, a re- 
cent experiment in two-dimensional systems with cavities 
showed clearer signals for the configuration with electric 
and magnetic fields perpendicular to each other, while in 
the parallel one the observed resonance feature presented 
puzzles lii 

Systems of particular importance are heterostructures 
forming a two-dimensional electron gas (2DEG), such as 
GaAs semiconductors. Here, the Rashba SOI is linear in 
momentum and provides an effective internal magnetic 



field about which the spin precesses. Realistic 2DEGs, 
moreover, contain disorder leading to momentum scat- 
tering which is responsible not only for the finite electric 
resistance but also for spin relaxation due to randomiza- 
tion of the internal fieldii22*22, Unlike for the conventional 
paramagnetic settingi^, no microscopic study of EDSR 
in such systems is available, but would be highly desir- 
able, also in view of the recent experimental activities. 
However, the interplay between SOI and disorder in or- 
bital space can be quite subtle, as was fully appreciated 
only recently in a number of experimental and theoretical 
studies on spin currentsii£*i2i2242ii22iS There, it turned 
out that the intrinsic spin Hall effecti^ in GaAs 2DEGs 
due to Rashba SOI does not survive the presence of dis- 
order, which, technically speaking, results from an un- 
usual cancellation of vertex corrections Conse- 
quently, spin currents in these systems are dominated by 
other (extrinsic) effects li^fiSi^ Similar concerns, for in- 
stance, apply to conclusions reached for EDSR in clean 
systemSiiS Thus, the role of disorder needs to be exam- 
ined carefully, and, in doing so, we show here that EDSR 
survives impurity scattering but acquires a line shape 
that is determined by disorder and SOI. 

We consider a geometry as shown in Fig. ^ where the 
electric rf field and the static external magnetic field are 
parallel and both in the plane of the 2DEG. At resonance, 
the spin polarization (magnetization) acquires a non-zero 
out-of-plane component. This we show first for a clean 
2DEG by deriving an effective Bloch equation. Turn- 
ing then to 2DEGs with disorder we treat the electric rf 
field in linear response and obtain for the magnetization 
a Lorentzian resonance whose width is given by a gener- 
alized D'yakonov-Perel spin relaxation rate. In addition, 
we find a shift of the resonance due to disorder and SOI 
which gives rise to an effective g-factor that depends on 
the magnetic field. Using a standard diagrammatic ap- 
proach to treat SOI and disorder systematically, we find 
that it is the spin vertex correction (coming from disor- 
der) that leads to the resonance, in stark contrast to zero 
frequency where the spin vertex vanishes^ Assuming re- 
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FIG. 1; Setup for electric dipole induced spin reso- 
nance (EDSR). The disordered two-dimensional electron 
gas (2DEG) is placed in the xy-plane (green area). An a.c. 
bias Vbias(t) with frequency uj generates an electric a.c. field 
E(t) along Gy which is in-plane and parallel to the externally 
applied magnetic field bo = g^s^o/'^ (assumed to be static). 
As a consequence of the orbital motion of the electron induced 
by E(f), the Rashba spin orbit interaction generates an 'inter- 
nal' magnetic a.c. field bi which is in-plane and perpendicu- 
lar to the momentum of the electron and thus to bo (ignoring 
disorder) . This configuration of fields is now analogous to con- 
ventional paramagnetic resonance and allows manipulation of 
the electron spin. 

alistic system parameters we identify the most promising 
regime for the experimental observation of EDSR. Fi- 
nally, we discuss the spin Hall current and show its rela- 
tion to EDSR. 



I. THE MODEL 

The 2DEG consists of non- interacting electrons of mass 
m and charge e which are subject to a random impurity 
potential V. In addition, we allow for a general SOI 
r2(p) • <T linear in momentum p and a static external 
magnetic field Bq applied in-plane, as well as a time- 
dependent electric field E(i) applied as a bias along Bq, 
see Fig. ^ The Hamiltonian for this system reads 

^ " i (P- ^A)' + fl(p).<T+(bo + bi(t)).rT + T/, (1) 

where hi{t) = — f {-^{t)) is the rf part of the internal 
'magnetic' field induced by the electric field via SOI, and 
A{t) — — c /* dt'E{t') is the associated vector potential, 
c being the speed of light. The Zeeman term contains 
6o = 5^bBo/2 and the Pauli matrices cr*, i — 1,2,3, 
while ri(p) is the 'zero-field spin splitting' due to the 
internal SOI field. 



II. EFFECTIVE BLOCK EQUATION FOR THE 
CLEAN SYSTEM 

We show now that for the described setup the dy- 
namics of the electron spin exhibits resonant behaviour 
(EDSR) generated by the electric rf field E(i). Starting 
with the simple case of no disorder {V — 0), we derive 
a Bloch equation for the spin dynamics (for weak SOI), 
from which the EDSR property immediately follows. We 
begin by noting that the density matrix p{t) is diagonal 



in momentum space, and its elements can be expanded in 
the spin basis as p(p,t) = J2t=o^'^(P'^)^'^ ^ where cr° — 
1. The expectation value of the spin is then {a^{t)) — 
Jd^p s'(p,t)/(27r)2 = S^{t), where Js°d^p/{2Tr)^ = 1 
due to normalization. [Henceforth, we refer to 5*(t) as 
polarization, becoming the magnetization when multi- 
plied with the Bohr magneton fiB-] In momentum space 
the symmetry is broken by the small SOI term such that 
the coefficients decompose into an isotropic and a small 
anisotropic part, s*(p, t) = s'^{p,t) + As^{p,t). Averaging 
the von Neumann equation for p(p, t) over directions of 
p we obtain 

§(p,i) = |(bo + bi(<))xs(p,i), (2) 

where we have dropped the angular average of ^^(p) x 
{s{p,t) + As(p,t)) since it is higher order in the SOI. 
Eq. ^ is now recognized as a Bloch equation describ- 
ing spin resonance. Indeed, specializing henceforth to 
Rashba SOI fi(p) = ap x e^,, where is a unit vector 
along the confinement axis, and taking E in the plane 
along bo||ej, results in the standard resonance setup 3& 
with an oscillating internal field bi(t) (x E x _L Bq. 
Tuning E(i) to resonance, i.e. to frequency lu = lul, 
with lul — q^bBq being the Larmor frequency, the spin 
starts to precess around the x-axis (in the yz-plane) with 
a Rabi frequency ojr = bi/Ti = eEoa/htUL given by the 
amplitude of the electric field Eq . 

III. POLARIZATION OF THE DISORDERED 
2DEG IN LINEAR RESPONSE 

Having established the existence of EDSR for the clean 
system we turn now to the realistic case of a disordered 
2DEG. For this we assume a dilute random distribu- 
tion of sho rt-ranged scatterers with the disorder average 
V{x)V{x') = (mr)~^5(x — x') taken to be i5-correlated 
and proportional to the mean free time r between elas- 
tic scattering events. The interplay between SOI and 
disorder can then be characterized by the dimensionless 
parameter x = 2ppaT/h measuring the precession angle 
around the internal field 2ppa/gfiB between scatterings. 

For A = and V — the eigenenergies of H become 
Es{p) — p'^ /2m + s bcs{p), with s = ±1 and the effective 
magnetic field bofr(p) = ap x -f bo. In the correspond- 
ing retarded (advanced) Green functions 

) 2 £^^E-Es{p)±i/2T' ^ ' 

the disorder manifests itself as a finite self-energy term 
i/2r generated by the disorder average. In the self- 
consistent Born approximation r is independent of mo- 
mentum due to the short range nature of V. We have 
also checked that a renormalization of the Zeeman split- 
ting, i.e. matrix valued corrections, are of order bg/Ep 
and can be neglected as the Fermi energy Ep is taken 
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to be the largest energy scale. Thus, the averaged Green 
function is given by ||2J) with the standard isotropic and 
spin independent term ij^T?^ 

We turn now to the explicit calculation of the spin po- 
larization (magnetization/zi^) S'^(lo) (per unit area) at 
frequency induced by the Fourier transform of the elec- 
tric field E(a;). Working in the linear response regime we 
start from the Kubo formula for S'*(w) averaged over dis- 
order and evaluate it using standard diagrammatic tech- 
niques. In Refi2^ such a calculation was performed for 
the static (w = 0) and zero field (Bq = 0) case, where, as 
a simplifying feature, the spin vertex correction turned 
out to vanish. For finite frequencies, however, this is no 
longer the case, and, as we shall sketch now, it is this 
vertex correction which leads to a finite out-of-plane po- 
larization S'^(lS) at resonance w = w^. To be specific, for 
<^ Ep the polarization becomes 
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(4) 



where summation over repeated indices is implied. Here, 
the velocity operator Vj = (i/h)[H, Xj] = Pj/m + a{ez x 
cr)j contains the usual spin-dependent term, and Tr — > 
J cP q / {2tt)'^ ti denotes momentum integration and trac- 
ing over spin states. In Eq.l^J we introduced the spin 
vertex correction S(a;) determined by the diagrammatic 
equation given in Fig. [3 




FIG. 2: Diagrammatic equation for the spin vertex 
correction S(tj). Iterating the equation for S*(tj) generates 
the ladder diagrams accounting for the dominant contribu- 
tion of the impurity scattering (denoted by a cross) to the 
Kubo formula Eq. I^J. The impurities causing momentum 
scattering of the electron are assumed to be short-ranged, 
isotropic, spin-independent, and randomly distributed in the 
2DEG. Here, ct* denotes the Pauli matrix for the ith compo- 
nent of the electron spin, G^^^{E) is the retarded (advanced) 
Green function averaged over disorder (see Eq. O, and r is 
the mean-free-time between scattering events. 

where the cross denotes the insertion of a factor 1/mr. 
The class of diagrams generated by E corresponds to the 
ladder approximation with an accuracy of order 1/ppl, 
I = ppT/m being the mean free path length. Thus, weak 
localization corrections j^iS being higher order in l/ppl, 
are not considered here. 

To further evaluate Eq.Q) we calculate S in the (de- 
composed) form S*(ti-') = X]i/=o ^*'^(^)''^'' ^'^^ ^^'^ limit- 
ing case of a strong magnetic field, bo ^ app, and for 
Bq = 0, respectively. In both cases(cf. the methods 
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FIG. 3: EDSR and linewidth in a disordered 2DEG. 

The reduced out-of-plane polarization M'^ — x{l — 1/A)E'^^ 
(cf. Eqs.Q, ®) is plotted as a function of the Larmor fre- 
quency ujl ~ gi-iBBo, with resonance peak at = a; of width 
F given in Eq. 0. Modulus (full line), real (dashed) and 
imaginary part (dotted) are shown for lot = 1 and x — 0.15. 
Inset: Peak height of at resonance as a function of lot for 
X = 0.2 and x = 0.1. 



section), the polarization is given by 



S^{uj) = eEj{uj)av2T€^j^ 



A 



, (5) 



where tiju is the Levi-Civita tensor, \{uj) = 1 — iwr, 
and V = denotes the two-dimensional density 

of states. The form of S(a;), however, is to be modi- 
fied accordingly due to the presence (or absence) of the 
magnetic field Bq. 

First, we consider the simpler case without magnetic 
field, i.e. Bq — 0. Here, the off-diagonal elements I]'^''(w) 
vanish. Thus, since the spin vertex is diagonal, we see 
that according to Eq. © S(tj) lies in-plane and perpen- 
dicular to E(a;), i.e. S'^{lo) vanishes identically for all lu 
in this case. In the considered geometry, with the elec- 
tric field along the y-axis, we are thus left with a finite 
in-plane polarization S^{lo) in the x-direction only. Fur- 
ther, for w — > the factor 1 — I/A vanishes and thus the 
spin vertex Y?^^ (0) drops out from Eq. |SJ . This way, the 
result oS^ is recovered. 

Next we consider the opposite case of a magnetic field 
Bq that is large compared to the internal SOI field. We 
characterize this regime by the small expansion param- 
eter a = app/2bo — x/2ujlt ^ 1, being the ratio of 
precession angles around the internal and external mag- 
netic fields, respectively, between scatterings. The field 
Bq leads to an equilibrium polarization along the y- axis 
due to Pauli paramagnetism given by Scq — vg^BB()/2. 
As we shall see, the EDSR response generated by the 
electric field E(i) will be proportional to S'eq, and thus a 
finite equilibrium magnetization is crucial. [We note that 
the Kubo formula, Eq. Q), describes the deviation from 
equilibrium only in response to the electric (and not the 
magnetic) field.] 
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With both, the internal and external magnetic fields 
present, the dispersion relation -Es(p) is no longer 
isotropic as 6eff depends on the direction of p. This com- 
plicates the momentum integrations considerably, and as 
a further complication, S becomes off-diagonal through 
the additional spin terms in the Green functions. For- 
tunately, for a <^ 1 analytical progress is still possible, 
and can be evaluated explicitly as we outline in the 
methods section. 

The polarization S'^{lu) is again given by Eq. 
but now with the spin vertex found for Bq (see 
Eq. HA3|I below), which has the finite off-diagonal ele- 
ments S^'^ = — E^^. Accordingly, the out-of-plane com- 
ponent is proportional to the reduced magnetization 
= xT,^^{l - 1/A) which is plotted in Fig. |3|as a func- 
tion of the Larmor frequency lul , exhibiting a pronounced 
peak at the resonance lul — u (up to a small shift, see be- 
low). This clearly shows that the EDSR resonance found 
in the clean case survives weak disorder, and, technically 
speaking, stems from the spin vertex correction. 

IV. EDSR AND LINEWIDTH 



nal and one by the internal field. At resonance and for 
ui^T ^ 1 the dynamical narrowing effect due to the exter- 
nal field is complete, but not the one due to the internal 
field, since in the linear response regime considered here 
the Rabi frequency is restricted to the regime ujrt <^ 1 
(see appendix). Thus, in order to describe complete dy- 
namical narrowing in the presence of disorder one needs 
to go beyond linear response, which, however, is outside 
the scope of the present work. For sufficiently strong E- 
fields in the limit cort 3> 1 it is plausible to expect that 
disorder becomes irrelevant for the spin dynamics. In this 
case, the Bloch equation (jJl applies and, like in standard 
paramagnetic resonance, the linewidth will then be given 
by the internal field bi (see also SecO. 

Next, the bare location ojl of the resonance gets renor- 
malized by the disorder leading to a small shift Suj <C lol, 
determined by the real part of 7 in Eq. ljA3p . 



Suj = - ^Re 7(w = Wi) = — — -. (8) 



Let us now analyze the line shape of the peaks close 
to resonance lul — in more detail, i.e. for \lul + 
w| <C \ujl - uj\ (left) and \u)l - <C \ujl + w| (right 
peak), respectively. In this case and to order Eq. (0) 
can be rewritten as a sum of two Lorentzians, with the 
transverse components (z = 1,3) becoming 



1 — iujT 



1 



1 



lul — uj + Slo — iT ujl + lu — Slu + iV J ' 



(6) 



where we denote loi = ujl for the in-plane {i = 1) and 
Lu^ = —iLu for the out-of-plane component (« = 3), respec- 
tively. Analogously to conventional paramagnetic reso- 
nance, the prefactor of can be rewritten in terms of 
the Rabi frequency uir = eEoa/hajL, i-e., the amplitude 
of the 'internal' rf field bi(t), and the equilibrium polar- 
ization — due to Pauli paramagnetism. 

Eq. 10 is the main result of this work. Here, F is the 
linewidth of the resonance peak and is explicitly given 
by the imaginary part of the damping function Eq. ljA4p 
evaluated at the (bare) resonance. 



1 



2ujlt'' 



■ Im 7(w = lol) = 7^-1 
zr 



2[l + (c^Lr)2 



(7) 

Here, the prefactor x'^ /2t is the familiar D'yakonov-Perel 
spin relaxation rat o^^i^? coming from internal random 
fields induced by disorder. The second factor accounts 
for the magnetic field. We note that F saturates for 
ujlt ^ 1 at the finite value x^/2r, and, thus, dynamical 
narrowingSi^ due to fast spin precession is incomplete. 
This reflects the SU(2) nature of the spin with two pre- 
cession axes on the Bloch sphere, one given by the exter- 



We note that for ujlt ^ 1 we get F ~ 56u}. Also, 
the shift 6u! depends non-monotonically on the magnetic 
field (wi), giving rise to an effective g-factor, gefj = 
g(l-ha;2/4(l + (ujLTf)), which is B-field dependent. Ob- 
servation of this effect in real measurements would pro- 
vide useful additional evidence for EDSR. Finally, the 
longitudinal component (i.e. along Bq) of the polariza- 
tion vanishes identically, i.e., S'^{uj) ~ 0. This is not quite 
unexpected in the linear response approximation, since 
it is known from conventional paramagnetic resonance^^ 
that changes in the occupation probability are nonlinear 
in the driving field. 

Let us now give some numbers to illustrate the EDSR 
effects discussed here (for details see Sec. |d1i. To quan- 
tify the amount of spins resonantly excited, we introduce 
the relative polarization P = (N^ — iV|)/(iV| -I- N^) given 
by the ratio of >5''^(a;) and the electron sheet density. Here 
and denote the average number of spins pointing 
up and down along the z— axis, respectively. Choosing 
now parameter values typical for a GaAs 2DEG sample 
of a few hundred micrometers in lateral size, a magnetic 
field of one Tesla, and an a.c. bias of 0.1 Volt and 8 
GHz we get P — 10~^ at resonance. This corresponds to 
Ni — Ni = 200 excess spins in the laser spot of a typ- 
ical optical measurement scheme,^^ and although being 
small, this is within reach of detection. Furthermore, the 
associated Rabi frequency, relaxation rate, and resonance 
shift, respectively, become a'i^/27r = 12 MHz, F/27r = 50 
MHz, and 6uj/2-k = 10 MHz. We note that the polariza- 
tion P increases at resonance with increasing frequency 
ijj (cf. inset in Fig. Eland Eqs. Q and iP). Thus, the 
higher the resonance frequency (and the magnetic field) 
the larger the EDSR signal. 
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V. SPIN HALL CURRENT 

We finally turn to a brief discussion of the associated 
spin Hall current for a homogeneous infinite 2DEG (for 
details see Sec. 10). Using the Heisenberg equation of 
motion we can express the spin Hall current, defined as 

— ({'7^, Wa;})/2r^ in terms of the transverse polariza- 
tion components as 

I'i^) = ^ [tujS\uj)+ujLS\co)] . (9) 

The importance of this relation lies in the fact that it 
establishes a direct link between the spin Hall current 
and the spin polarization, and that the latter quantity is 
experimentally accessible by known measurement tech- 
niques. We note further that for if = Eq. would 
describe circular motion of and in the x-z plane 
(at resonance), i.e. 'pure' precession of the spin around 
the y-axis. A finite spin current, therefore, character- 
izes the deviation from pure precession and gives rise to 
nutation of the spin. From the ratio hlf/2Ey the as- 
sociated a.c. spin Hall conductivity follows, which, at 
resonance, is explicitly given by (see Eq. cr^jj'^'* — 

{e/4:Tr)iujLT/{l + 2A). Thus, for the high field regime 
ujlt ^ 1 we see that the spin Hall conductivity assumes 
a universal value a^'y^'^ = |e|/87r. Quite remarkably, this 
is the same value as obtained in the absence of magnetic 
fields and at finite frequencies 

In conclusion, we have investigated the spin polariza- 
tion of a 2DEG in the presence of spin-orbit interaction 
and disorder. We have shown that carrier spins in 
a specific field configuration with an electric rf field 
give rise to EDSR with a linewidth coming from spin 
relaxation due to disorder and spin-orbit interaction. 
Our results emphasize the importance of tunable SOI 
for coherent spin manipulation by electric means in 
semiconductors . 
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APPENDIX A: METHODS 

Here we sketch some of the important steps of the cal- 
culation. A more detailed account can be foimd in Sec. 

ini 

Calculation of the spin vertex correction. For the 

evaluation of Eq. we introduce the spin-spin (X^") 
and the spin- momentum (Y^^) terms, resulting from the 



spin and the normal part of Vj, respectively, 
X^-'iLu) = J-TT{G^{q,Ep) G''iq,EF + hu;) a''] 

y« (c.) = -^Tr f G^(<z, Ep) G^(q, Ep + ^ | 

(Al) 

and expand the spin vertex in terms of Pauli matrices, 
= Y^^^qT?" {uj)a'' with the coefficients S'"' be- 
ing represented by a 4 x 4-matrix. Here, the Fourier 
transform of a function /(t) is defined as /(w) = 
J^^ dtf{t)e'-'^y2'K. With this notation we then obtain 
from Eq. Q for the polarization 

nh 

(A2) 

We proceed using the identity (5"" S^^ = 
'^l'^Y^f_L=Q^ai3^^'a' such that the diagrammatic equation 
in Fig. [21 factorizes and can be solved. This way, the 
matrix elements of the "diffuson" are obtained in the 
form W{uj) = [[1- Xiujyf^Y" and the polarization 
S^{ijj) is expressed in terms of the 4 x 4-matrices F, X, 
and S, which were evaluated for the two cases bp = 
and 6o 3> app, respectively. 

The case Bq = 0. Performing the integrals in 
Eq. for Bo = we find X^°{uj) = 1/A and 

X^'^{u}) — A/(A^ -I- x^), while the in-plane components 
are equal, = = {X^° + X^^)/2 (these results 
have already been obtained inS in the evaluation of the 
spin Hall conductivity). The off-diagonal components 
vanish so that the spin vertex follows simply as 
— (1— X''^)^^. The spin- momentum matrix is found 
as y« = -ae^j-a/A + 0{a^) (see also Sec. El). 

The case Bq 7^ 0. In the case of the strong magnetic 
field, i.e., Bq ^ it turns out that, in leading order in the 
SOI-strength a the spin-momentum term keeps the same 
form as in the absence of the magnetic field, i.e. Y^^^ = 
— ae^ys/A. The polarization is thus again described by 
Eq. (0), however, with the spin vertex modified by the 
presence of the field Bo (cf. Eq. (|A3I) ). In this case, the 
calculation of S(ti') is more involved since, starting from 
weak SOI-strength a ^ 1, the leading order power in a 
changes precisely at resonance. Thus, an expansion up 
to second order in a of the spin-spin term X with the 
subsequent matrix inversion E = [1 — X)^^ is required 
to consistently describe the shape of 5* in terms of the 
spin vertex at resonance. 

As a result (see Sec. IbJ, we obtain 

S](c.) = 

I ih \ 

^ , ' 

(A3) 
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where the imaginary part of the complex function 



7(w) 



3(A - 1)A3 - (a;^r)^A(6A - 1) - {loltY 
2A(A2 + (c^Lr)2)2 



(A4) 



characterizes the hnewidth induced by the impurity scat- 
tering via the SOI, as described in the main text. The 
functions x^^ii(ui) and x^7i3(u;) (exphcitly given in Eq. 
(|B11|) ') are neghgible for obtaining the polarization S'^{ijj). 
However, for obtaining the spin Hall current they are 
needed since the the leading terms of T}^ and Y?"^ cancel 
each other in Eq. ljE5|l . 



APPENDIX B: CALCULATION OF THE SPIN 
VERTEX CORRECTION 

Here, we give a more detailed account of the evalua- 
tion of the polarization S^{lo)^ Eq. (jA2|) . in terms of the 
spin-spin and the spin-momentum terms, X, Y , respec- 
tively, and the spin vertex correction E. The simultane- 
ous presence of the internal and external fields, a(p x e^) 
and bo, resp., breaks the symmetry in orbital space such 
that no closed analytical expression for the integrals in 
Eq. (|A1I) and hence for X^^'^ can be obtained in the gen- 
eral case. However, the most important regime for EDSR 
is given by the regime where the internal field is much 
smaller than the (perpendicular) static external magnetic 
field, as in standard paramagnetic resonance (see also 
next section). Thus, without any essential restriction we 
can concentrate on the regime with the SOI being small 
compared to bg, i.e. a — app/'^bo — x/2lolt <^ 1. First, 
upon inspection of Eq. (|A2p we note that the contribu- 
tion of Y to the polarization is due to the momentum- 
part of the velocity and thus must vanish in the absence 
of SOI. Thus, the leading order term in Eq. (jA2|) coming 
from Y is at least linear in a (assuming analyticity). More 
precisely, with a calculation similar to the one outlined 
below for X, we obtain for the spin- momentum diagram 



W3 



A(c.; 



-f 0(^3), 



(Bl) 



i.e. the same result as before for bo — 0. Note that 
only odd powers in a appear due to the symmetry con- 
straints induced by the angular integration occurring in 
Y. Thus, the expression Eq. lO is linear in the SOI 
a in leading order. In order to expand the polarization 
Eq. in leading order in a cx a it is therefore suffi- 
cient to calculate the spin-spin diagram with setting a 
to zero and retaining only the bo dependence. This way, 
we obtain the spin vertex correction S(cl;) which is singu- 
lar at resonance, i.e. when lu — lul (Larmor frequency), 
refiecting the presence of Rabi oscillations. This shows, 
however, that at resonance the next-to-leading order con- 
tributions of the spin-spin diagram X = X(o) + a^X(^2) 
become relevant^ for the matrix inversion and must be 
kept. Indeed, they represent the dominant contribution 



in S(w) = {l — X(^a){uj) — a'^X(^2){^)) ^ if the determinant 
of the first term 1 — -^(o) vanishes. Obviously, at reso- 
nance the dominant a-dependence of S becomes 
Hence, we concentrate on the evaluation of the spin-spin 
diagram up to order of a^, with A^(2) characterizing the 
behavior of the polarization around the resonance (where 
our analysis is valid). 

The spin-spin diagram is given by 

1 



(l+W-Q2_s5(Q)+ir/2) 

1 

(1 - 02 _ s'B(Q) - ir/2) 



expressed in terms of the dimensionless quantities^ w = 
TiLij/Ep, r = h/EpT, Q = q/pp, Bq — bo/Ep and the 
effective magnetic field 



B(Q) 



boff(QpF) 



Ep 



Bn 



2aBo(Q X e^) (B3) 



with modulus 



B(Q) =Bo\/l + 4aBo- (Q x e,)+Aa?Q-^. (B4) 

Here, Bo = Bo/i?o is the unit vector of the external field 
taken along the y-axis such that the mixed product in 
Ea. (jB4|) becomes Bo ■ (Q x e^) = — Qcos^s. 

Trace. In Ea. ljB2p we introduced the trace over spin 
states 

TZ'{^) = tr{'T^(l + s B(Q) • a)a''ll + s' B(Q) • a)} 

= 4<5^''[,5^°<5,,,, +J^^°<5,,_,,] 

+ 4 [(5''^"<5''"S^(Q) + (5'^o^'^^"B,(Q)]s5,,,, 

+ 4ss'B^(Q)S,(Q)(5'^^"(5''^o , (B5) 

where 5'^^'^ = 1 — (5^° etc., and where summation over 
repeated indices is implied. There are terms containing 
none, one, or two normalized magnetic fields B(Q) = 
B(Q)/_B(Q), which is relevant for the momentum inte- 
gration. The trace T"^^, (Q) and the direction of Bq de- 
termine the matrix structure of the spin-spin diagram, 
i.e. which components X^^ are nonzero. 

Momentum integration. The components X^^'^ are 
obtained by the momentum integration in Eq. ljB2|) 
where the poles Qi of the denominator yield the main 
contribution. Assuming that Ep represents the largest 
energy scale such that w, r, and are small compared 
to one, the poles are located essentially at \Q\ ~ 1 (cf. the 
denominator in Eq. ljB2p 'l with corrections 0{'W, BQ,r). 
Making use of a ^ 1 we expand the denominator of the 
integrand in Eq. I|B2|) in a up to second order. From this 
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ImQ 



+1 



""Q-i ReQ 



FIG. 4: View of the poles Qi of the retarded and advanced 
Green functions in the complex plane (cf. Eq. ijBBfl ). The 
contour Ci of the momentum integral (Eq. IBSH running 
from to +00 is shown. 



we obtain the following four poles 



Qi = 


+1- 


h ki 


r 


O2 = 


-1 - 


Vk2 


r 

— I — 
2 


Qi = 


+1- 


^k3 


r 

-'2 


Q4 = 


-1 - 


F ki 


r 

+ ^2 



(B6) 



where 



w Bq 2 . 2 

Ki^2 = i-^ T s-^ — h -Bqis cos 93 =F a s_Bo sin 

^3,4 = "Fs'-^ + Boas' cos 93 =F a^s Bq sin^ . (B7) 

These poles are approximately of order one with a small 
correction and a small imaginary part r /2 — h/2EpT, 
thus showing that the above approximation for small a 
is self-consistent. Decomposing into linear factors the 
spin-spin diagram can be recast into the form 



dQQT^Z'(Q) 



(B8) 



where {■■■)^ denotes integration over the polar angle (nor- 
malized by 27r), and c^p = 1 + o^Bq^s + s') sin ip. 

Note that the Q-dependence of \/B{Q) in T^^^, for 
a generates additional singularities at Q^/f, = 
exp (±i(/5)/2a when analytically continued into the com- 
plex plane away from the real axis. The application of 
complex contour integration is thus non-trivial. A direct 
calculation, however, carried out by decomposing the de- 
nominator in Ea. (|B8|) into partial fractions shows that 
these poles do not contribute to X within the accuracy 
0{w,B^,r). 

The subsequent angular integrations become simple 
when expanding in the small parameter a. In this way. 



we obtain the components X^^" and via the matrix in- 
version S = (1 — X)~^ the spin vertex correction in Eq. 
l|A3p of the paper. In particular, the vertex components 



_ {ujltY + (A - l)A + a;27n 



with the complex damping function 



(B9) 



characterizing the linewidth and the functions 



7ii(Lj,a;L) = - 
7i3(w,wl) 



2ujltX^ 

[(wlt)2 + A2]2 



(Bll) 



are relevant for the subsequent calculation of the spin po- 
larization and spin Hall current. The frequency depen- 
dence and resonance behaviour of the spin polarization 
and current are discussed in the main text. 



APPENDIX C: REGIME OF VALIDITY 

We now give a summary of the parameters controlling 
the regime of validity of the present theory. 

A first constraint ensures the validity of the linear re- 
sponse approach. For this, we give a heuristic argument 
based on the analogy to conventional ESR^ expressed 
by Eq. ||2Jl of the paper. For this case, we consider the 
Bloch equations 



^i=7[Sx (Bo + Bi(t))], 
^3=7[Sx (Bo + Bi(i))] 



5^=7[Sx (Bo + Bi(t))]2- 



3 T2 

S^-S, 



eq 



Ti 



(CI) 



describing magnetic moments subject to a constant mag- 
netic field Bo 1 1 62 and a circularly polarized field Bi(i) 
perpendicular to it, which oscillates with frequency uo. 
The familiar steady state solutions in the rotating frame 
for the longitudinal component S"^ and the transverse 
components S*" and S"" , respectively are 



5* — Sq. 



1 + Aw22.2 



5" = S 



UJ1T2 



cq 



(C2) 
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where uji = 7B1, lul = jBq, and Auj = ojl — i^- 

The resulting transverse polarization close to res- 
onance is, thus, proportional to loiI[{u} — lol)"^ + 
Lo\{Ti/T2) + I/T2] with the phenomenological relaxation 
rate 1/T2. Thus, two relaxation terms are present, viz. 
the 'external' damping given by l/T-f and an intrinsic 
term w^(Ti/r2) given by the driving rf field itself. 

Similarly, the same intrinsic mechanism should be ex- 
pected if the driving field Bi is generated by a SOT 
mediated bias like in the case considered in the paper. 
We thus anticipat e a total spin relaxation rate of the 
form^ \/uj\ + with Rabi frequency lor ~ eEoa/huJL 
derived at resonance from Eq. (0). Here, Eq denotes the 
amplitude of the electric field E(i) = EoSy cos(wi) and F 
is given by Eq. Q. However, the Rabi frequency occur- 
ring in the rate \/ujji + F^, being E-field dependent, must 
be negligible for a polarization which is calculated in 
linear response with respect to E(i). This imposes the 
self-consistent condition 



UeE, 







< 2apF 1 -I- 



1 



2\l 



(C3) 

for the validity of the linear response approach. A more 
systematic approach for estimating the validity of the 
linear response regime requires an explicit evaluation of 
the non-linear response, which, however, is beyond the 
scope of the present work. 

Secondly, in order to carry out the momentum integrals 
in Ea. ljB2|l we introduced a condition limiting the SOI 
strength 



hujL 



< 1. 



(C4) 



This constraint not only simplified our analysis but also 
defines the most interesting regime for EDSR. Indeed, in 
order to have a pronounced resonance, the width of the 
resonance peak needs to be smaller than the resonance 
frequency, i.e. F <C ti^L, which is equivalent to appx ^ 
hojL (see Eq. {7}). For self-consistency we need to assume 
X < 1 (see the text before Eq. jnj), and thus we see that 
a ^ 1 ensures F <C wl- 

In this context we note the somewhat counterintuitive 
fact that the height of the resonance decreases with in- 
creasing SOI, see Eq. jnjl. Indeed, on one hand the po- 
larization is proportional to a via the driving rf field, and 
thus increases with increasing SOI. On the other hand, at 
resonance the polarization becomes proportional to 1/F 
(due to disorder) which gives then rise to a suppression 
factor Thus, in total the polarization decrease as 

1/a with increasing SOI at resonance. 

Our last constraints 



bo huj h 
Ep ' Ep ' EpT 



< 1 



(C5) 



correspond to the physically relevant situation where the 
Fermi energy Ep is the largest energy in the system. Fur- 
ther, the condition x = 2appT/h <ti 1 does not restrict 



Description 


Parameter 


Value 


sheet density 




4 X 10^^ cm~^ 


effective mass 




0.067 me 


scattering time 


T 


2 X 10^" s 


frequency 




8 GHz 


Larmor frequency 


f = fr = uJr /2-K 


8 GHz 


R,ashba Parameter 




10~^^ eVcm 


electric field 


E 


1.66 Vcm~^ 


polarization 


P 


10"* 


SOI vs. scattering 


X 


0.1 


spin relaxation rate 


r/2-K 


0.05 GHz 


resonance shift 


5lo/2-k 


0.01 GHz 


Rabi frequency 




0.012 GHz 


validity conditions 


linear response 




0.27 


relative SOI strength 




0.05 



TABLE I: Numerical estimates 



the validity of Eqs. ljA3|l and (|A4p but permits us to rep- 
resent Eq. (IBJ in terms of two Lorentzians. In the case 
LOT K, 1, however, it becomes equivalent to the inequality 
llC4ll . 



APPENDIX D: NUMERICAL ESTIMATES 

To illustrate the predicted effects we now evaluate the 
polarization explicitly using typical GaAs parameters (cf. 
table P) , thereby making sure that we stay within the 
range of validity of our approximations. With a typical 
sheet density n2 = p|,/27r?l^ = 4 x 10^^ m~^, effective 
mass m* — 0.067 me and a high mobility scattering time 
T = 2 X 10"^^ s taken fromiS we can estimate the maxi- 
mum polarization P as the ratio of the peak polarization 
per unit area and the sheet density 



93 

p max 

n2 



eEm* 



{2Trhn2)^2aT 



1 



2(l+t^V2) 



(Dl) 



In order to stay within the condition (jC4|) we choose a 
small Rashba - parameter a = lO"^"* eVm and find x = 
0.1. Assuming a realistic microwave frequency w = 50 x 
10^ s^^ corresponding to = 1 and a voltage amplitude 
of y = 0.1 V over a sample length of / = 600 fim we find 
an electric field Eq = 166 Vm^^ and a polarization of 



P = 10~ 



(D2) 



9 



Note that the size of the chosen i?— field satisfies the lin- 
ear response condition HC3|I (and poses no severe limita- 
tion for a real experiment). 

The corresponding number of excess spins -/V| — Ni in a 
laser spot of size 5/im x 5/im is 200. This number is mea- 
surable with state-of-the-art optical detection techniques 
such as Faraday rotation.^* 

We can further quantify the peak width T and the 
frequency shift Scu. Making use of Eqs. {TJ and (jHJ of 
the main text we find 



r = 0.3 X 10^ s"^ 
Sio = 0.06 X 10^ s"^ 



(D3) 



As a further characterization of the resonance we esti- 
mate the Rabi frequency lob., given by the amplitude of 
bi{t) in Eq. ((SJ. Assuming a bias E(t) = E^ey cos{u!t) 
we find 



^ £|2^ ^ 0.08 X 10« s-\ 



(D4) 



evaluated at resonance ■ cj^ = ui with the parameters 
given above. A summary of the above calculation and a 
check of the constraints Eqs. (|C4IC3|I is given in table 



APPENDIX E: SPIN HALL CURRENT AND 
POLARIZATION 

We show now that the obtained polarization 
(magnetization/zis) 5** can be related to the spin current 
(defined below) via an exact relation. More generally, we 
consider the spin density operator 



pXx) = ^K,<5(x-X)}, 



(El) 



defined as the (symmetrized) product of the spin with 
the particle density operator (5(x — x) where x is the 
position operator. Integrating over space (homogeneous 
limit) and taking expectation values we get the spin 
polarization S** = / cPx < p*(x) >. The spin cur- 
rent density associated with is defined in the usual 



waj* 



13.16.30.31.32.33.41 



jl(x,i) = -K,j,(x)} 



(E2) 



in terms of the current operator jfc(x, t) = ^{S{x—ic),Vk} 
where, in contrast to the linear response treatment of 
the paper, the velocity operator Vk = i/h[H,Xk] = (pk — 
(e/c)Ak)/m + a{(T x e^);; contains the kinetic momentum 
including the (homogenous) vector potential A. 

The two operators and jl are related via the Heisen- 
berg equation of motion 



(E3) 



given by the Hamiltonian Eq. . Analogous toi^ where 
the Rashba- and Dresselhaus SOI has been considered it 
forms an exact operator identity 



2am 



(x, <) + V • j 1 (x, t) - - (x, t) 



dt 



- [p2(x,t)feo,2 -P^(x,<)6o,j,] 



-p" (x, <) + V • (x, t) = —jI (x, t) 

2 

- - [p^(x,t)&o,^-p\x,t)6o,z] 

^p=^(x,<)+V-j3(x,t) = [ji(x,t)-fj2(x,i)] 

2 

- ^ [p\x,t)&o,a - p^(x,t)&o,a.] , 

(E4) 

for the case of an additional static magnetic field with 
components 69. i — x,y,z, which holds independently 
of the impurity potential as commutes with the posi- 
tion operator. 

In deriving Eq. (|E4|I the definition of arises natu- 
rally as a divergence term of a current associated with 
the spin density. Together with the time derivative it 
forms the left-hand side of a continuity equation. The 
right hand side, however, is nonzero and describes the 
dynamics of the spin due to the external magnetic field 
bo and the internal SOI field. The definition of Ea. (|E2l) 
as a 'spin current' is thus ambiguous'^'''^^ and it is not 
clear to what extent the quantity Eq. (|E2|) can be iden- 
tified with actual spin transport, i.e. with spin polarized 
currents which are experimentally accessible.*'^ 

In spite of the above concerns we note that in the ho- 
mogeneous limit the spin Hall current can be expressed 
entirely in terms of the polarization. Namely, going over 
to the spin Hall current 11 = J < jl{x) > such that 
the gradient in Eq. (|E4|) vanishes we find the expectation 
value of the spin Hall current given by 



A, 



(E5) 



[This relation can be obtained directly from the Heisen- 
berg equation of motion da^ /dt = i[H, cr^], and by noting 
that /3 = (j^Px ■] Since S'^ vanishes for = it is obvious 
from Eq. (jE5|) that there is no spin Hall current in the 
dc limit w — !■ for a homogenous infinite sample. 
This means a generalization of the argument given inli^ 
to the case of a finite magnetic field. For finite frequen- 
cies, however, Eq. (jESp predicts a non- vanishing oscil- 
lating spin current expressed in terms of the polarization 
components perpendicular to the applied electric rf field. 
With the results for inserted we find the ac spin Hall 
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conductivity evaluated at resonance (w = ujl) as 



^3, res 



iujlt 



Ey 



(E6) 



We emphasize that this relation provides a direct link 
between the experimentally accessible polarization and 
the spin current. For w^r ^ 1 this becomes 



giving a universal value for the spin current at resonance. 
It is quite remarkable that the same result, Eq. (|E7|I . 
can be obtained when inserting the solutions S'* obtained 
from the Bloch equations Eq. IjCip close at resonance 
into Eq. (|E5)| . 



3, res 
'xy 



e 

'8^ 



(E7) 
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As a general property of linear SOI the first order in a 
vanishes due to the symmetry in the angular integration. 
Indeed, we note that the angular dependence (in the in- 
tegrals in Eq. IIAH ') comes from the magnetic field (Eq. 
lB3l l) where ip always occurs in terms of a trigonometric 
function simultaneously with a factor a. Expanded in a the 
linear terms thus vanish upon angular integration. 
*5 In this Section]^ the capital letters B,Bo and B,Bq de- 
note dimensionless magnetic fields measured in units of the 
Fermi energy Ep. 
assuming Ti = T2 for simplicity 

corresponding to a magnetic field B ~ 1 T for \g\ = 0.44 



